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ABSTRACT 


Natural gas, a non renewable source of energy can be found trapped inside the earth’s surface in form of gas hydrates. Recovery of gas from class 1 gas 
hydrate reservoirs using depressurisation techniques can be simulated using the reservoir-tank model. This paper extends the analysis of the Oregon site 
using this model. It studies the effect of deliverability constant on the gas recovery. The paper also argues that shutting down of the well at non equilibrium 
pressures actually leads to pressure build up in the well due to hydrate dissociation which in turn increases net recovery of gas from the well. We analyse 
sequential closing and opening of the well at different dissociation rates for the hydrate as a method to maximise the net recovery of gas from the well ina 
fixed time span. 
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1. INTRODUCTION 


Non renewable sources of energy have been in continual demand for over three centuries now. Natural gas or Methane is a 
widely used non renewable source of energy. Natural gas and oil are generally found together below the earth's surface 
where the gas is trapped between the cracks of rocks. Another form in which natural gas occurs is in the form of hydrates. A 
natural gas hydrate is a methane molecule trapped inside water molecule cage. The abundance of such reserves is almost 
four times more than the known coal reserve. Large abundance of hydrate reserves and higher calorific value of methane are 
going to play an instrumental role in making natural gas hydrates next big thing in Energy industry. The hydrates are generally 
classified into four classes on the basis of how they occur inside the surface. Class 1 hydrates are those hydrates where a 
layer of hydrates exists above a free natural gas layer. In the present study we have tried to study recovery of gas from a 
class 1 hydrate and how to increase the recovery from such a hydrate. 


2. MODEL 


The model was applied to a site in Oregon which was in the category of class 1 hydrates. A potentially significant source of 
methane hydrate lies eighty kilometres off the coast of Oregon, USA (Trehu and Fleuh, 2001). The site is evaluated using the 
tank reservoir model (Khatanair, 2002). The value for different parameters useful in the modelling and their values for the site 
are mentioned in the Table 1. A homogeneous and isotropic hydrate layer is assumed to be on a free gas layer. The reservoir 
is assumed similar to a tank where the gas obtained from well is comprised of initial free gas as well as gas obtained from 
dissociation. The dissociation is considered analogous to ice melting-it occurs along the front inside of the entire volume 
(Ahmadi 2004). The tank reservoir model assumes that the reservoir behaves as a tank and there is only one centrally 
located production well. This reservoir is assumed to be homogenous and isotropic with a hydrate zone and a free gas zone 
beneath. The driving force for dissociation is only depressurisation. In the hydrate zone there are three phases assumed - 
hydrate (solid), natural gas (gas) and water (liquid). Instantaneous equilibrium in terms of pressure and temperature is 
assumed throughout the gas zone. The reservoir behaves as a closed system with no mass or energy transfer through its 
boundaries. Lastly it is assumed that there is no free water in the reservoir and all the water present is produced due to 
hydrate dissociation. In a typical class 1 hydrate there is a layer of hydrate over a layer of free gas. In the tank model, let the 
hydrate layer be of height hug there is only single hydrate phase present in this layer and it is assumed that all the water and 
free gas produced goes to the lower layer. At any instant the height of hydrate layer is hy changed from the initial hydrate 
layer length of hxi. In a typical class 1 hydrate there is a layer of hydrate over a layer of free gas. In the tank model, let the 
hydrate layer be of height hig there is only single hydrate phase present in this layer and it is assumed that all the water and 
free gas produced goes to the lower layer. At any instant the height of hydrate layer is hg changed from the initial hydrate 
layer length of hrgi .The hydrate decomposes at the equilibrium pressure of Peg , which at the isothermal conditions of 6°C is 
5400 KPa and the water and gas influx from this dissociation is We and Gen. The molar rate of dissociation of hydrate is (Kim 
et al, 1987). 
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Table 1 
Nomenclature 


A Reservoir area (m?) 1500 
Bg Gas formation volume factor . 
Bgi Initial gas formation volume factor : 
Bgh Hydrate formation volume factor * 
Bw Water formation volume factor - 
Gp Gas produced (m°) - 
Gen Gas influx from hydrates(m°) - 
Gii Initial free gas in place(m’) - 
Grr Remaining free gas (m°) - 
Gui Initial hydrated gas in place (m°) - 
Gur Hydrated gas remaining(m®) - 
hy Hydrate zone thickness(m) - 
hui Initial hydrate thickness (m) 115 
hgi Initial gas zone thickness (m) 300 
Pi Initial reservoir pressure, (KPa) 5400 
Pwt Flowing wellbore pressure, (KPa) 1600 
Og Gas production rate, (m/s) - 
qw Water production rate, (m%/s) - 
Sw Water saturation in the gas zone at timet = 
Swi Initial water saturation 0.15 
Wop Water produced, (m°) - 
Wen Water influx from hydrates, (m°) - 
o Porosity 0.2 


Where, K, = 2.2838 L0** scf gas/acre-psi-day 
And == 16950°R = 9416.666K 


Also, The water influx is given by WeH=0.00083KGen 


The free gas layer contains free gas and water and the saturation is Swi initially. As the hydrate decomposes, this layer 
increases to hg and saturation becomes Sw. The gas at any point in the system from dissociation of hydrate is Gen. The 
pressure at any point in this layer is P. The gas and water is continuously removed from the layer through a well. The well 
flowing pressure is Pwr and is 1600 KPa. The water produced from the well is Wp and the gas produced is Gp. 


The rate at which gas is obtained from the well is modelled as 
_ dtp 


_ 
@={F=cx |e —F..<° (Cis in m%/s/KPa) 


The pressure P of the free gas is slowly decreased to this pressure Pw. 


The rate at which water is obtained from the well is 


= > _ Betas 
i - pele a 


The whole tank is assumed to have a surface area of A and porosity of the equations is reached by the understanding that in 
the whole process the volume of the system is conserved. The changes in volume of the three phases are considered 
separately and the net change in volume is equated to 0. 


Change in the volume of Hydrate = AVH= (Gui-G-H1)BH 
= Ab(1-SwiA hy 


Change in the volume of Water = A Vw=(Wp-Wen)Bw 
= Ad [Swix hg -Swx (hg+A hg)] = Ad(Swi- Sw)hg 


wwered 


Change in the volume of Gas = Volume change due to pressure conditions change + Volume change due to gas production 
Volume change due to pressure conditions change = Gri (Bgi-Bg) 

Volume change due to gas production = (Gp-Gen)Bg 

Total Change in volume of Gas AVg=Gri(Bgi-Bg)+(Gp-Gex) Bg 

Total volume change of Reservoir = 0 


disc@™very 
Corer the uncovered 


AVu+ A Vw +A Vg = 0 


> — (Gui}-Gui})BuH + (Wp-Wen)Bw+ Gri(Bgi-Bg)+(Gp-Gen) Bg = 0 


3. METHODOLOGY 


The volume changes of hydrate, free gas and water amount give the total change in volume of the system. In the reservoir 
model this change of volume is assumed to be zero and thus one equation is obtained. Differentiating this equation with 
respect to time gives an equation explaining the variation of Pressure. Further, the volume formation factors were calculated 
using gas laws and the kinetics equation which was taken from Kim et al (1987). The system thus could be modelled using 
three coupled differential equations in P, hg and Swi. These equations were then solved in MATLAB. The initial conditions for 
solving are taken as follows 
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Table 2 


P= 5400 KPa , hg=300 m ,Swi= 0.15. 


Pressure at the end of every 100 seconds in an open well 


(C=0.05 m%/s/KPa) 


Results of the simulation are discussed in the next section. 


on 4. SHUT DOWN 


4474 The results of the section above show that the pressure in the well falls very rapidly with 


4020 the time for which the well is opened. The reservoir pressure at the end of first 100 


3563 seconds of opening the well is found to be 4934KPa. The fall in the pressure at the end 


3146 of every 100 seconds when compared to initial pressure is reported in the Table 2. If a 


2735 well at a pressure less than the equilibrium pressure is currently not used to obtain gas, 


2356 it is said to be in shut down. This shutting down of the well can actually increase the net 


2024 recoverable gas from the well by promoting decomposition of the gas hydrate above the 


1767 free gas layer. Pressure difference from equilibrium pressure is the driving force of this 


Table 3 


1621 dissociation. The decomposition of the hydrate continues until the equilibrium pressure 
is re-established in the well, which is 5400Kpa for given conditions. The equations for 
pressure and volume variation during a shut down were solved in MATLAB with initial 


Optimal set Ci for which recoverable gas is maximised (Ci in conditions being the conditions at the end of first 100 seconds for an open well. The 


m*/s/KPa) 


variation of pressure with the time passed was studied and plotted. It takes around 12000 
days for the pressure to come very close to the equilibrium pressure, without ever 
100 actually becoming equal. It is practically impossible to shut the well down for 12000 days 


200 and then reopen it again thus continuing the cycle. However, the gain in total recovery 


300 due to shut down for a short time is also considerable and can be utilised when 


400 maximising the overall recovery in a given time span. 


10500 5. INCREASING RECOVERY BY SEQUENTIAL OPENING 


eno AND CLOSING OF THE WELL 


20700 The results in the previous sections show that pressure is build up in the well if it is shut 
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20800 down after opening it for some time. The initial equilibrium pressure is retained. For a 

typical well, the pressure at the end of first 100 seconds is 4934KPa .If the well is shut 
down thereafter; the time to regain the equilibrium pressure is around 12000 days .The shutting down though doesn't 
contribute to gas flow but builds the pressure which increases the recovery if the well is subsequently opened. Though this 
sequential opening and closing of the well is practically infeasible, we should not under appreciate the increase in the gas 
recovery due to a shut-down. The effect of shut down and corresponding pressure build was studied when the gas obtained 
from the well for a given time is optimised. Let’s say we have to make 10 choices of C for a given well. If C is non zero that 
means well is opened for 100 seconds with C equal to the chosen value. On the other hand we can choose C to be zero and 
assume the well is shut down for 10000 seconds and the rate of dissociation is 10° times more than normal rate. The rate can 
be increased by either decreasing the energy barrier (E) by catalysis or increasing the temperature (T) of the well or both. So 
we can either choose to operate the well for 100 seconds at a non-zero value of C or shut it down for a 1000 seconds and a 
increased dissociation rate by a factor of 10°. 


(0,0.1] 100 Open 
0 10000 Shut down 


In other words, we need to choose Ci for i=1 to 10 .The selection should be such that the total gas recovery from the well is 
maximised. The code is implemented in MATLAB where a random initial guess of Ci for i=1 to 10 is generated using randi(). 
This guess is then improved such that the obtained gas from the well increases with every improvement in the guess. Finally 
an optimum set of Ci for i=1 to 10 is returned by the obtained gas maximizer function. To obtain a more reliable optimum set 
of Ci the process is repeated a large number of times. The optimum set of Cj that is obtained maximum number of times is 
reported as the final optimum set (Table 3). 


6. RESULTS 

1. The Graph 1 shows the fall in the reservoir pressure with the time for which the well has been left open. The initial 
pressure is equilibrium pressure at 6*0C which is 5400 KPa.It can be seen that only within 1000 seconds which is only 
16.67 minutes ,the reservoir pressure reaches the well pressure of 1600 KPa 


2. The well deliverability constant(C) also has an effect on the drop in the pressure of the reservoir. The Graph 2 shows that 
for larger values of C the pressure falls more rapidly, as is expected. 


3. The pressure in a closed well as a function of the time (in days) is shown in Graph 3. The equilibrium is re-established 
after only around 12000 days which is a very long time when compared to 100 seconds for the pressure to fall to this 
value. 


4. The initial guess for Ci=0.01 randi (0,10). The guess is improved continuously by the fval function in MATLAB so that the 
obtained gas is maximised. 


5. The Ci values for which the total gas is maximised are given in Graph 4. 


7. CONCLUSION 


The analysis explains a lot about the nature of the class 1 hydrates and also help to look into future ways to increase the gas 
recovery from a given site. The study shows that it is almost practically infeasible to obtain significant continuous flow of gas 
from a class 1 hydrate at the given site .The shutting down of well increases the pressure build up inside the well. Therefore, 
the total recovery from the well can be increased by sequential opening and closing of the well where during the shut-down 
hydrate dissociates and increases the total free gas in the well. The dissociation rate during shut down should be increased 
appropriately to make the recovery significant thus making the process economically feasible. 
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Figure 1 Figure 2 
Reservoir pressure versus the time Effect of C(in m*/s/KPa) on the pressure in the reservoir 
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Figure 3 Figure 4 
Pressure built up in a closed well Sequential closing and opening of well to optimise gas recovery.(Ciin 
m*/s/KPa) 
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